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Abstract. We study a mathematical model of a blood circulation network which is a 
generalization of the coronary model proposed by Smith, Pullan and Hunter. We prove the 
existence and uniqueness of the solution to the initial-boundary value problem and discuss 
the continuity of dependence of the solution and its derivatives on initial, boundary and 
forcing functions and their derivatives. 

1 Introduction 



In a recent paper ||13|| , Smith, Pullan and Hunter propose a mathematical model of blood 
circulation in the coronary network, and conduct a numerical analysis. In their model, major 
vessels (with cross-sectional areas larger than a certain value) are treated as a connected one- 
dimensional network, and small vessels, such as arterioles, capillaries and venules, are treated 
as lumped elements which are connected to the network of vessels. The flow on vessels are 
assumed to be incompressible, homogeneous, Newtonian, and has a small Reynolds number. 
Thus, the mass balance equation and Navier-Stokes equation can be written to describe the 
pressure and the flow rate on vessels. Equations on lumped elements are written in analogy 
with the current and voltage in an electric circuit. The result is an initial-boundary value 
problem of a system of hyperbolic type partial differential equations coupled at junctions 
of the network. Although the result of the numerical analysis conducted in |1| matches 
closely with measured data, the well-posedness problem of the system of partial differential 
equations, that is, the existence, uniqueness and the continuous dependence on initial and 
boundary data of the solution, has not been established before. The main objective of this 
paper is to establish the well-posedness. We prove that the system is well-posed under 
certain natural conditions. This work is an extension of our earlier work [jEJ on a model 
of blood circulation in the brain. The main differences between the two models are that 



the network configuration in |y| is more complicated owing to the presence of Willis loops, 



but the coupling junction conditions in the model of |U| are more complicated due to the 
different formulation and the inclusion of the capillaries and veinal system. We combine 
both features in a more general system with the hope that our result will be useful in the 
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modelling of circulation systems of higher complexity, including the whole body circulation 
system. 

Before stating our system, let us briefly describe the model in ||13|| . Let Pi and Ri 
represent the pressure and radius on the i-th vessel, respectively, and let Vi be the cross- 
sectional average of the axial component Vi^ x of the velocity on the i-th vessel. Assuming 
that the radial component t>; jr of the velocity is small compared to the axial component u ia . 
of the velocity, one can write equations of mass balance 



dRi T dRi RidVi 
+ V~ + 



and momentum balance 

dVi 



dt 
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Here v is the viscosity constant and 
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is the energy quantity. Taking into consideration of no-slip boundary condition (vi )X = 
if r = Ri), the viscous axisymmetry (dvi, x /dr = if r = 0), and the fact that V{ is the 
cross-sectional average of Vi jX , Smith, Pullan and Hunter propose the velocity profile 



Vi-. r, x 



7^ + 2 

li 



Vi(x) 



Ri 



ii 



where 7 i is a positive number. Using this profile and the mass balance condition ( |1 . 1|) , Eq. 
(II. 2ft becomes 
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2vai Vi 



IR 



1.3) 



with 



7, + 2 

li 



e (l,oo). 



The pressure P, and the radius i?j are related by a function 

P i = P i {x,R i ). 



In [13 , it is assumed that 



Pi (x, 



C 



Ri, 
Ro 
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where C, Rq and (3 are constants. We do not make such an assumption, only assume that it 
is a different iable function and 

dPi 

w, >0 

for all x and Let 

At = 7i Rl Qi = AiVi 

be the cross-section area and the flow rate, respectively. It can be shown that the system of 
equations ( |1.1|) and Ql.3| ) is equivalent to 

dA i + dQ 1 Q 
dt dx 

+ on— — r + 



<9t dx \ A J p dx cti — I Ai 

Also, by rescaling the spatial variable x, we may assume that each vessel is parameterized 

to x e (0,1). 

The system of differential equations are supplemented with the initial condition 

Pi(x,0) = Pi (x), Q i (x,0) = Q I l (x) (1.5) 

and boundary conditions. Boundary conditions at each end of the vessel are given according 
to the type of the end. If it is an external end of the network, either the pressure 

Pi = Pi (t) (1.6) 

or the flow rate 

Qi = Q?(t) (1.7) 

is specified. If the end is a branching junction, a junction connecting several vessels, let 
ji, . . . ,j u and j u+ i, . . . , jfj, denote the incoming and outgoing vessels, respectively. One im- 
poses the mass balance condition 

V /J, 

1=1 l'=u+l 

and the momentum balance condition 

Pjl = Aj l (Pj t — Pjunc) ) Pj v ~~Qf = Aji> {Pjunc ~ Pj ir ) (1-9) 

for I = 1 . . . , v, I' — u+1, • • • , /i, where p i are small positive constants and Pj Unc is the pressure 
at the junction. (In [fC^ , every branching junction connects only three vessels, it is either a 
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bifurcation point of one artery into two smaller ones or a joining point of two veins into a 
bigger one. Our prescription allows more general configuration of the network, including the 
presence of Willis loops.) If the end is a transitional junction, which connects the vessel to a 
network of arterioles, capillaries and venules, we follow the so called microcirculation model 
proposed in ]13|, [yj. Generalizing from |jl3| , arterioles or venules connected to the vessel ji 



are represented by a lumped resistive element Rj r The capillary bed is also represented by 
a resistive element Rc- Rj^s are connected to Rc through capacitive elements C\ and C 2 
on the two ends. 




Figure 1: Microcirculation model of the network of arterioles, capillaries and venules. 

Let ji, . . . , j u be the arteries and let j u+ i, . . . , be the veins that are connected to a arteriole- 
capillary- venule network. The boundary conditions for P JP Qj l are 



RjiQji (1, t) = P n (1, t) - P Cl (t), for I = 1, . . . , v, 
Rj v Qj v (0, t) = P C2 (0, t) - P jv (t) , for V = v + 1, . . . , fi 



and 



HP HP 

1=1 l'=u+l 



dt ^— ' ' dt 

where Pc,; i = 1, 2 represent the pressure in the capacitive elements C\, C 2 , and 



;i.io) 



Qc= PCl p ^ (1-12) 

represents the flow rate in the resistive element Rc- (In [13], there is only one artery and 
one vein connected to the system of arteriole-capillary-venule at the two ends. We do not 
rule out the possibility of multiple arteries and veins join together to such a system.) 
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The system we study in this paper consists of the equations 



m, 0Q i _ r 

-+- (Jj Ji, 

° OX 16 (0,1), t >0 (1.13) 

dt 1 dx C% dx l ' 



and the initial and boundary conditions given by ( |1.5|) - (|1.12|) . For convenience, we also use 
the vector form 

(Ui) t + Bi {Ui) x = Fi (1.14) 
where Ui = (Pi,Qi), F { = (fi,gi) and 



Bi 



di 
bi la 



Eq. ( |1.4| ) is a special case of this system where 

dP t A { aCfi fdPA- 1 aQ t aQ\ dA t AnisaQ, 



1 1 p A} \dAj ' 1 A t ' ' ^ ^ 9x a-lA' 

We do not assume any particular form of these functions though, they are general differ- 
entiable functions of (x,t,Pi,Qi). Our basic assumptions are a, > and Aj > e for some 
positive constant 6q. Other assumptions will follow. Apart from the junction conditions, 
this system is the same as the one we study in [12]. Also, the junction conditions in |12| is 
the special case of ([L.8|) -( |L9|) above with p i = 0. As in |12|], we use a fixed point principle 
to prove the solvability of the problem. Substituting a pair of functions (pi, qi) for (Pj, Qi) in 
the coefficients a iy b iy Ci, Ai and forcing functions fi, g iy the system becomes linear. That is, 
all the functions a,, etc. are independent of unknowns. If the linear system has a unique so- 
lution, then, one can establish a mapping from (p^, to the linear problem solution (Pj, Qi). 
If one also shows that this mapping has a unique fixed point, then the fixed point is neces- 
sarily the unique solution of the quasilinear system. Hence, we shall first give a condition 
for the linear system to have a unique solution, then examine under what conditions the 
mapping has a unique fixed point. The first aspect of the problem is investigated in Section 
2 and the second in Section 3. We also prove a result on the continuity of dependence of 
solutions on the initial, boundary and forcing functions for linear and quasilinear systems, 
thus, completing the analysis of the well-posedness of the problem. In spite of similarity in 
parts of the analysis to the one used in ||12|| , the more general branching junction condition 
and the new transitional junction conditions require more careful treatments. Hence, there 
are substantial variations in the analysis. For completeness and to benefit the reader, we 
include all the major arguments in this paper. 
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2 The linear system 



In this section, we analyze ( |1.13|) as a linear system with a«, bi, Ci, fi, gi, Ai independent of 
Pi and Q{. The initial and boundary conditions are given by (lL5|)- (P^D except that the 
junction condition ( |1.9|) is substituted by the more general condition 



Pji Q£ ~ Aj t (Pj t Pjunc) + Cj t , Pj v — (Pj unc Pj ir ) + Cj v (2.1) 

where Ci are differentiable functions of (x, t) . The inclusion of Ci is needed in the next 
section in order that the result of this section can be extended to the quasilinear system. 
We give conditions for the linear system to have a unique global solution. The conditions 
are most naturally given in terms of the eigenvalues of the matrix Bi, which have the form 

Af = ^ + Ui, Af = ^ - Ui, 

where 

Ui = \Jcj + Oib~. 

These eigenvalues are real if 

c- + CIA > 0, x G (0, 1) , t > 0, i = 1, . . . , n. (2.2) 

In this case, 

Af(M)>0, Af (x,t) < Af (x,t) (2.3) 



and the system is hyperbolic. Under the condition ( p.2|) , we show that the linear system has 
a unique solution if 

Af (0,t) < 0, Af (l,t) < 0, i = l,...,n 

which is equivalent to 

> 0, t > 0, i = l,...,n. (2.4) 
at x = 0, 1 only. It needs not hold for x G (0, 1). 

Theorem 2.1 Assume that the functions a i} bi, Ci, fi, gi, Ai and Ci are independent of 
{Pi-iQi}- Suppose that these functions and the initial and boundary functions Pf , Q\, P Z B , 
Qf all have bounded first-order derivatives. Suppose also that ai > 0, Ai > and that the 
conditions (\2.Q) and (\2. A) hold. Then, for any T > there is a unique solution in a bounded 



subset of the space C([0, 1] x [0,T],R ) to the linear system ( \1.13{ ) with the initial and 



boundary conditions given by ( fQ} j-flO | j, and ( \2.1\ ). 
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Proof. We first show that the system has a unique solution for < t < 5 for some 5 > 0. 
The proof is based on the method of characteristics and a fixed point principle. For systems 
defined on only one branch with boundary conditions of the forms of ( |1.6| ) or ( |1 . 7| ) , this is a 
standard approach. In our case, special care is needed to handle the junction conditions. 

Consider the z-th branch. From any point (£, r) on the left, right, and lower boundary of 
the rectangle D =: [0,1] x [0,T], we construct the left-going and right-going characteristic 
curves x — xf (£; £, r) and x = xf (i; £, r) by 

Af (xf,t), xf (t)=£, 
Af (*f,t), xf (r)=e, 

respectively, where Af and Af are the two eigenvalues of the matrix Bi. By the uniqueness of 
solutions to these differential equations, a left-going (resp. right-going) characteristic curve 
cannot intersect with another left-going (resp. right-going) characteristic curve. Let Xf and 
Xf be the right-most left-going and left-most right-going characteristic curves, 

x = xf (t; 1, 0) and x = xf (t; 0, 0) 

starting from the lower boundary of -D, respectively. It can be shown from (|2.3| ) that the 
two curves can have at most one intersection. Let £j be the value of t at the intersection. 
If the two curves do not intersect in D, we simply define tj = T. By condition ( |2.4j ), X/' 
cannot reach the right vertical line x — 1 at any £ > 0, and by Af > 0, Xf cannot reach the 
vertical line x = at any £ > 0. Thus, the rectangle Di =: [0, 1] x [0, ti] can be divided into 
three parts 

where Df is between the vertical line x = and the characteristic curve Xf , Df is between 
the two characteristic curves, and Df is between Xf and x = 1. 



dxf 
dxf 




Figure 2: Three parts of Di 
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We show that there is a < t{ such that the solution (Pj, Qf) for the i-th branch exists in 
the restriction of Di to the strip {0 < t < Si}. 

First, observe that the initial conditions alone determine the solution completely in the 
central region Df . This follows from the theory of first-order linear hyperbolic systems and 
the fact that from any point (x,t) G Df, the two characteristic curves, followed backwards, 
must land on the horizontal line t = 0. (The latter is a consequence of (|2.3|) .) To extend 
the solution to other parts of Di, we make a change of unknowns and derive a set of integral 
equations. Note that if =: (— Aj , af) and if =: (— Af , af) are the left eigenvectors of Bi 
corresponding to Af and Xf , respectively Introduce new unknowns 

r, = ifUi = -Af P t + OiQi, Sl = ifUi = -Af P + a 4 Q 4 . (2.5) 

The system ( |1.13|) can be written in terms of and Sj by multiplying the left eigenvectors 
to ( p,.14|) and substituting in 

P = ^- (n - 8i ), Q< = 2^ - A ^ s • ( 2 - 6 ) 

This results in the equations 

dfn = F* (x, t, n, Si ) , dfsi = Ff (x, t, n, Si ) , (2.7) 

where 
and 

if (x, t, r i; Si ) = /f P + (9f If ) U t , Ff (x, t, r u Si ) = ifFi + (dflf) Ui. (2.9) 

(A differential operator acting on a vector means that it acts on each component of the 
vector.) Let (x,t) G P«. We integrate the first equation of ( p.7|) along the right-going 
characteristic curve x R (t; r) which passes through (x, t) and reaches the left or lower 
boundary of Di at (£, t). It can be shown that for (x,t) G Df U Df, r = 0, and for 
(x,t) G Df, £ = 0. In the former case, we obtain 

r t (x, t) = r[ (0 + / Pf (xf (f; C, 0) , t>, n, Si ) df (2.10) 
Jo 

In the latter case, we have 

n(x,t) =n(0,T)+ [ F t R (xf{t';0,T),t',r l ,s l )dt'. (2.11) 
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Similarly, by integrating the second equation of (|2.7|) along the left-going characteristic curve 
x i (t] r ) that passes through both (x, t) and (£, r) (which is on either the right or lower 
boundary of Di), the equations are 



Si (x, t) = s{ (0 + / if (xf (*'; f , 0) , t', n, Si) dt' 
Jo 



if (x, t) e Df U Df and 



(x, t) = Si (1, r) + / if (xf it'- 1, r) , t', n, s<) dt' 



(2.12) 



(2.13) 



if (x, i) G .Df 1 . These are the integral equations we need. 

For any Si < ti we use -D^., Df s . and Z?^ to denote the restrictions of Df, Df and Df 
to the strip {0 < t < Si}, respectively. First, consider the case where the end of the branch is 
an external end. We discuss the case of a left end only, the case of a right end can be treated 
similarly. If the boundary condition is given by ( |1.6| ), we define Sj = Sj/s where e < 1 is any 
constant. Using the first equation of Q2.6|) in the integral equations (|2.11|) and (|2.12|) , 



/ 



Ti {x, t) 
Si (x, t) 



2 Ui (0, r) Pf (r) + es ( (0, r) + J if (xf (*'; 0, r) , t', n, esj dt 



-*[(£) + - / if (a£ (t';e,0),t',r i ,es i )dt' 



) 



(2.14) 

This is a fixed point equation for (r», s$) if we define the right hand side as a mapping of an 
operator K on (r»,Si) in a bounded subset of C (Df s _ U Df s .,M. 2> ). In a standard approach, 
it can be shown that K is a contraction mapping if Si is sufficiently small. Hence, the fixed 
point exists and is unique, and the solution (r i; sf) is uniquely extended to Df s U Df s .. If 
the boundary condition is given by (|1.7|) , we define Sj = Si/e, where e > is so small such 
that 

< 1, r G (U, . 



Af (0, r) 



The fixed point equation is then 
/ 2aiUi (0, r) 



n 0, t) 
Si (x,t) 



V 



) 



Af(0,r) -y ' ' Af(0,r) 
-4 (£) + -/ Ft (4 (f; C, 0) , t', r», esi) dt' 



(2.15) 

By a similar argument, the solution can again be uniquely extended. 

We next extend the solution to either Df g or Df s if the end is a branching junction. 
In this case, we shall extend the solution on all the branches that are connected to the 
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same junction simultaneously. Let ji, . . . ,j u be the incoming and j v +i, ■ ■ ■ , j^, the outgoing 
branches to the junction. Equations ( |1.8|) , ( |2.1|) and ( |2.6| ) give rise to a 2/i x [i homogenous 
system of linear (ordinary) differential equations for (1, t), Si(l,t), i = ji,...,j v and 
Ti (0,t), Si (0,t), i = j v +i, . . . 



Pn 
An 


d Q 


(M) 


_ £li 

Ah 




(M) = 




c 

-£-Pi(l,t), i = j 2 ,...,j u , 


Ph 
An 


d Q 


(M) 


_ £h_ 

Ah 


~Pn 


(i,*) = 




C- 

- -j- - Pi (0, t) , i = j u+1 , ...J^ 


V 

















i=i 



l'=v+l 



(2.16) 

Differentiate the last equation with respect to t and regard Sj ± (1, t),..., Sj v (l,t), Tj v+1 (0, £),..., 
Tj (0, t) as unknowns. The derivatives of unknowns can be solved from (|2.16|) because the 



coefficient matrix of dsi/dt and dri/dt in Q2.16D , 



/ 



P.h X h ( 1 ' t ) 



2uj 1 a,j 1 Aj 1 

o 



\ 



2u jl a jl A jl (l,t) 



y 2uj 1 aj 1 (i,t) 



2u 32 a i2( 1 ' t ) 



2« jM % M (0,t) / 



has the determinant 



11 n,.n. A. h A 11 



Aa 

\ u n a 3i A n (M) u jV a frAfr (°'*) ^ ft, ' 

Since Af < < Af at the junction, the determinant is not zero. Thus, the derivatives of 
the unknowns, Sj l (1, t) and Tj , (0, t), are each a linear combination of the functions (1, t), 
Sj l (1, t), Tj , (0, £), s_j , (0, t) together with the derivatives of Tj x (1, t) and Sj , (0, £), Z = 1, . . . , u, 
V — v + 1, . . . , [A. Integrating and using the initial condition determined by ( |1.5| ) and (J2.5|) , 
we can write 

s 4 (1, r) = s 4 (1, 0) + K (?) r n (1, r) + ^ m} ;; (r) s,-, (0, r) + / H t dt', (2.17) 
z=i j'=i/+i ^° 

for i = ji, . . ., j v and 

r J (0,r)=r,(0,0) + ^n} ; (r)r j! (l,r)+ ^ nj, (r) Sj , (0, t) + / H.dt', (2.18) 

!=1 J'=i/+1 ^° 
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for i = jv+x, . . . ,jfi, where m*, n l - are continuous functions and Hi are linear combinations of 
Tj t (1, t), Sj t (1, t), Tj , (0, t), Sj , (0, t) and Cj (£) with coefficients depending only on t. Choose 
an e > such that 



£ max 



( r )l K ( r )| f < X > i = h,---,3» re [0,ti 



i=i i=i 



and introduce 

r • = — s • = Z = 1 i//' = z/ + l u 

Then, from (|2.10|) - (|2.13|) , the integral equations for the 2/i unknowns fj v sj t , jv,, §j,, I = 
1, . . . , v, I' = v + 1, . . . , // constitute a fixed point equation, u> = Kw, where 

w = if jn ■ ■ ■ > ?j v , Sj 1 , . . . , Sj v , Tj v+1 , . . . , 7j M , Sj„ +1 , • • • , % M ) (2-19) 

and 

Sjl (1, 0) + e (ELi m} k f jk (1, r) + ££ =i/+1 m)J jkl (0, r)) + £ H h dt> + £ Fj< . . . , 
r 3v+1 (0, 0) + e (ELi (1, r) + J2£>=u + i n)J 3y (1, r)) + J Q T if il/+ A + £ . . . , 

(2.20) 

in which 

F R = F R (x R t' er- s ) F L = F L (x L t' er ■ s ) 
for / = 1, . . . , v, and 

for I' = v + 1, . . . , //. It can be shown by a standard argument that K is a contraction 
mapping in the space 

V /J. 

x j= .\[c (//;.„ u Df u8i , r 2 ) x J] '* u n';^ , i 

Z=l l=i/+l 

if Sj is sufficiently small. Hence, it has a unique fixed point in Xj. This extends the solution 
(ri, Si) for the neighboring branches of the junction. 

It remains to extend the solution to a region adjacent to a transitional junction. Similar to 
the case of a branching junction, we simultaneousely treat all the branches that are connected 
to the same transitional junction. Let ji, . . . ,j v be the arteries and j v+ ±, . . . , j M be the veins. 
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The condition connecting the vessels are given by fll.lUj) , ( |1 . 1 1| ) and ( 1.12| ). Differentiate the 
equations in ( |1.10|) with respect to t, the resulting equations together with is a linear 

system of the derivatives of the functions (1, t), Sj t (1, t), r, , (0, t), Sj , (0, t), I — 1, . . . , v, 
V — v + 1, . . . , fj,, and Pci (t), Pc 2 The coefficient matrix of dsj t (1, t) /dt, drj v (0, t) /dt, 
for I — 1, . . . , v, I' — v + 1, . . . , ii and P' c , P' C2 is 



D 1 B 1 
D 2 B 2 
I 2 



where 



D 1 




R h X h + 1 



2u 



31 



+ 



2u 



'3" jv 



2u jv a ju 



+ 



2u 



3v 



R 3y. X t 1 



2u 3, a 3, 



2u j 

3p. 



I2 is the 2x2 identity matrix, and B\, B 2 are some constant matrices. Since all the elements 
of the diagonal matrices Di and D 2 are positive, the system can be uniquely solved for these 
derivatives. Thus, each of dsj t (1, t) /dt, drj v (0, t) /dt for I — 1, . . . , v, I' — v + 1, . . . , \i is a 
linear combination of (l,t), Sj t (l,t), rj, (0,t), Sj, (0,t), and dr^ (l,t) /dt, dsj , (0,t) /dt, 
I = 1, . . . , v, V = v + 1, . . . , /i as well as Pc x and Pc 2 - We can also eliminate Pq x and Pc 2 
by (|1.10| ). Integrating the resulting equations, we obtain equations (|2.17|) - (|2.18|) for some 
functions m*, n*, and Hi. The remaining of the previous paragraph can then be used here 
to give the extension of the solution to the left or right regions for the branches. 

Finally, if we let 5 be the minimum of all Si occurring above, we see that 5 > and the 
solution exists and is unique in (x,t) e D$ =: [0,1] x [0,5]. Observe that 6 depends only 
on the bounds of the system functions aj, etc., the initial and boundary functions P/, etc., 
and their first-order derivatives in D = [0, 1] x [0, T]. Hence, it is independent of t, and we 
can extend the solution successively in the time intervals [0,5], [5,25], etc. In this way, the 
solution is obtained in D in finitely many steps. ■ 

We next derive an estimate of the deviation of solution in term of the deviations of the 
initial, boundary and forcing functions. This estimate is needed in the next section. For 
any vector function v = (t>i, . . . ,Vk) defined in C (X; ~R k ), we use 

, where X represents a closed subset of either R or M. 2 



v \ x to denote the norm 



max,- 



\C(X) 



Lemma 2.1 Let U = (P,Q) and U = yP,Qj be two solutions of the linear problem ( \1.14 ) 
with different initial, boundary, and forcing functions. Suppose the conditions of Theorem 
hold for both solutions. Suppose also that there is a positive lower bound for all A; L . Then, 
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there exists a constant M > 0, independent of initial, boundary and forcing functions, such 
that 



u-u 


< M ( 


pi -pi 


+ 


Q 1 -Q 1 


+ 


pB _ pB 


+ 


Q B — Q B 




C(D S ) V 




C[0,1] 


C[0,1] 




C[0,5] 



C[0,S] 



+ 



C-C 



C[0,S] 



f-f 



C(D S ) 



+ 8\g-g\ 



C(D S ) 

(2.21) 



Proof. We need only prove ( [2.21j ) for a S < min^ {Si}, where Si represents the constants 
occurring in the proof of Theorem [2.1| . This is because for larger S, we can divide the 
interval [0,S] into subintervals, each has a length less than minj{£j}, and apply fl2.21|) in 
each subinterval. We can then take the maximum on each side of the inequalities to derive 
the inequality of in [0, S}. In the sequel, Df, Dg and Df are the restrictions of Df, Df and 
Df to the strip {0 < t < S}, respectively. 

By linearity, U — U is the solution of the system with the initial, boundary and forcing 
functions P} - P(, Qf - Q(, P* - Pf, Qf - Qf , h - k 9i ~ 9i ^dQ - C t . Let n, h, s h 
§i be defined as in the proof of Theorem |2.1| , corresponding to U — U . We show that these 
quantities have upper bounds in the form of the right hand side of (|2.21|) in Df , D L and 



In Df , (|2TT0| ) and (gig ) hold. Notice that the functions F* and are linear in r,-, and 
Si. Hence, there exists a constant M (we will use M generically for any constant bounds 
that are independent of solutions) such that 

Rf (t) + flf (t) < \ri\ cm + | S [| C[0)1] + mJ* {Rf it') + Sf (0 + Tf (i 1 )) df, 



where 



and 



Rf (t) = sup \ri(x,t)\, Sf (t) — sup \si(x,t)\ 

{x:(x,t)eDf} {x:(x,t)eDf} 



Tt (t) = sup 

{x:(x,t)€Df} 



fi (x, t) - fi (x, t) + \g { (x, t) - gi (x, t)\) . 
Hence, by Gronwall's inequality (see, e.g. M, p.327]), 



(2.22) 
(2.23) 



Rf + (t)<M[ K| C[M + K 



lc[o,i] 



S sup T t ° (t) 

te{0,8) 



for t G [0, S]. This proves that Rf and Sf have upper bounds in the form of the right side 
rf©- 
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We next consider the left or right regions if it is adjacent to an external end. Since the 
both cases are similar, we will only treat the case where the left end is an external. The 
integral equations to be used are ( |2.14j ) or ( |2.15| ) according to the type of the boundary 
condition. The resulting inequality has the form 



Rnt)+snt)<°snt)+M(\sf 



C[0,1] 



+ e 



B\ 



\C[0,S] 



+ / (R^{r) + Snr)+Tt{r))dr 



where £f is either Pf or Qf depending on the boundary condition, and Rf , Sf' and are 
defined in the same way as in fl2.22| )-( p.23| ), with Df substituted by Df U Df , and a > is 
a positive constant such that o = e if the boundary condition is (|1.6| ) and 



a = e sup 

te(o,<5) 



Af (0,t) 



Af(0,t) 



< 1 



if the boundary condition is ( |1 . 7| ) . Replacing M by (1 — a) M, we can write 



Rt (t) + S{< (t) < M (\sl 



C[0,1] 



IC[0,<5] 



+ 



RHt) + SHt)+THt) )dr 



Hence, by Gronwall's inequality 



Rf(t)+Snt)<M^\si\ cm + \^\ cm 



5 max T L (t) . 
te(o,s) ! 



This proves that both R\ (t) and (t) have upper bounds in the form of the right hand 
side of (P^TQ. 

We next extend the estimate to Df s or Df~ 5 if the end is either a branching junction or 
a transitional junction. In either case, the solutions on the branches j\, . . . , connecting to 
the junction constitute a fixed point of the operator K, which is defined in Q2.20Q . Let 

W(t) = j2 (R% (t) + S? (*)) + E (Ri (t) + % (t) 



i=i 



l'=u+l 



where Rf and Sj^ are defined as in ( gg|) with Df substituted by Df U Df. Then, from 
w = Kw and in view of the assumption that A, has a positive lower bound for all i and 
t > 0, we can deduce 



\l=l l'=u+l , 

+ M [E I r* I c[o,i] + E | s i | c[0)1] +J\w(r)+T (r) ) rfrj , 



, i=i 
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where 



1=1 



V=v+1 



1=1 



and (t) is defined as in (|2723|) with £>f substituted by Df U Df. Replacing M by 
(1 — cr) M, we obtain 



C[0,1] 



W(t)<M EKU,i] + EK / (W(r)+T(r))dr 

\l=l l'=v 

Hence, by Gronwall's inequality, 

W It) < M V Irf L rn +5 max T (t) ) . 



. «=1 



This leads to an upper bound in the form of the right hand side of ( |2.21 ) for Rf (t), (t), 
i = ji, ■ ■ and Rf (t), S[ (t), i = j v+1 , . . . , j M . 

We have thus obtained an upper bound in the form of the right hand side of (|2.21|) for 
the quantities |r$ — Ti\ c , D \ anc ^ I s * — s i\c{D s )- The conclusion of the lemma follows now from 



3 The quasilinear system 

In this section, we study the quasilinear system where the coefficients Oj, bi, q, fi, gi, A; and 
Ci depend on both (x,t) and (Pi,Qi). Under certain conditions, we show that the system 
has a unique local solution. We then present a theorem on the continuity of dependence of 
the solution on initial, boundary and forcing function. 

The basic idea in the proof of the existence of solution is to construct an iterative sequence. 
Substituting any vector function (pi,qi) for (Pj,Qi) in aj, etc., the system becomes linear. 
Thus, we can use Theorem |2J] to get a solution (Pi,Qi). This defines a mapping 5* from 
u =: (pi,qi) to U =: (Pi,Qi), and the solution for the quasilinear system is a fixed point of 
S. If there is a subset of a Banach space that is invariant under 5", then, we can construct a 
sequence 

u k+ i = Su k , k = 0, 1, — 

In the case where the limit exists and is unique, it gives rise to fixed point of S. This is our 
approach in this section. 

In this approach, conditions (|2.2j) and Q2.4p are repeatedly used. One might want to 



impose them for all the values of the variables. This would give the existence and uniqueness 
for the global solution, as in the case of the linear system. However, such a requirement is 
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so restrictive that even the original system ( |1.4j ) cannot meet it. Therefore, we will impose 
them only for t — 0, and obtain the local solution for the quasilinear system. 



Theorem 3.1 Assume that the initial and boundary functions P i , Q{, P t B , Qf and the 
system functions ai, hi, Ci, fi, gi, Ai and Ci all have continuous first- order derivatives with 
respect to each variable. Suppose that ai is positive and A% has a positive lower bound for all 
the values of their arguments, and that conditions ( \2.2[ )- ft2^\ ) hold at t = 0. Suppose also 
that the initial functions P[, Q\ satisfy any relevant boundary conditions at t = 0. Then, 
for some 5 > 0, there is a unique solution for < t < 5 to the quasilinear system ( \1.13j ) with 
the initial and boundary conditions given by ([!.$) - fil73() , (\2. 1\ ), and ( \1.10j) - ([i.lQ) . 



Proof. We first consider the simpler case where U 1 =: {P^Q 1 ) = 0. Let v = {^}, 
v i — (Pi, 1i) De a family of vector functions (not necessarily constitutes a solution) that 
satisfy the initial and boundary conditions. Substitute v for U in the functions a«, b iy Ci, fi, 



gi, Ai and C,. Then, the system becomes linear and we can invoke Theorem |2.1| to obtain 
a solution U to the linear system. This defines a mapping S : v t— > U. A solution to the 
quasilinear system is then a fixed point of S. We will choose a subset X$ t M of a Banach 
space such that (1) SXs,m C Xs,m , and (2) S is contracting in X$ } m - For any scalar 
or vector function / G C k (D$), let \f\ k $ denote the maximum norm of all the k-th order 

derivatives of / in D$. (If / is a vector function, \f\ kS = maxj l\fi\ k Let Cb (-D<5,K 2ri ) 
denote the subset of the vector-valued functions in C (Ds,M. 2n ) that satisfy the initial and 
boundary conditions. We seek Xs,m i n the form 

X SMo = {v e C B (D 5 ,M. 2n ) : \v\ 0>s < M , \v\ ljS < Mi| (3.1) 

where Mq is an arbitrary positive constant and M\ is a constant to be determined. Note 
that by the vanishing initial condition, for any M±, \U\ 1S < M\ implies |{/| s < M\5. Hence, 
for any M , we can ensure \U\ QS < M by reducing 5. It remains, therefore, only to show 
that for Mi sufficiently large and 5 sufficiently small, \vL s < Mi implies l^u^i < Mi. 
Throughout this proof, we use M to represent any positive constant that may depend on Mi 
but is otherwise independent of v and 5, and use M for any constant that is independent of 
Mi, v and S. The values of M or M in different occurrences need not be equal. 
Let U = Sv and let fj and Sj be defined by Q2.5Q . On each branch, we show that 

ms X {\(r i )J,\( Si )J,}<M 1 (3.2) 

and 

m ax{|(ri)J , \{si\\} < Mi (3.3) 
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in Df, D$ and Df if M\ is large and S is small. (Recall that Df etc. are the intersections 
Df fl -D«5 etc., respectively.) In fact, only ( |3.2|) needs to be shown. To see this, first ob- 
serve that the vanishing initial condition and the compatibility of the initial and boundary 
conditions gives 



max \ \P- 



lc[o,<5] 



,\Q 



B\ 



\C[0,S] 



< M5. 



Hence, we obtain from Lemma \2.1\ with U — that 

\UL S < M5. 



(3.4) 



From (57?) and (2J3), there are constants M and M such that 

\df ri \ < \l?Fi\ + \dflf \ \Ui\ <M + M5 } 



WsA < \li FA + |Sf Zf I |^| < M + M5 



(3.5) 



for each z = 1, . . . , n. Hence, Q3.3| ) follows from ( 
in We also note that ([T§ and (gTJ) imply 



( ^.5[ ) and the definition of df and 



d?UA n]t <M + M5 t \dfUA nx <M + MS 



(3.6) 



for all This will be used later. 

We first consider the middle region Df, where the solution (r,, Sj) satisfies the integral 
equations (|2.10| ) and ( |2.12|) with r\ — s\ = 0. Differentiating the equations with respect to 
x, we have 



(3.7) 



°t 

( s 4 = ffl.Ui (x,t) + [ [(l^F t ) x + {U i ) x - (lf) x (%U t )] (xf) x dt. 



Here, we used an identity from |5], p. 469]: 

d r b 

f(x(t),t)Dg(x(t),t)dt 

= f (x (b) , b) g x (x (6) , b) x ? (6) - / (x (a) , a) (x (a) , a) x 5 (a) 



(3.8) 



+ 



[f x (x (t) , t) Dg (x (t) ,t)-Df {x (t) , t) g x (x (t) , t)} (t) dt 



where x (t) is a function such that x (b) = £ and D = + x' (t) J^. (Notice that (6) = 1.) 
Let 

(3.9) 



i#(t)= sup {INJx,*)!}, = sup {\( Si ) x (x,t)\} 

{x:(x,t)eDf} {x:(p,t)eB°} 
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From 



and (|3.7|), we derive 

Rf (t) + Sf it) < M5 + M f (1 + Rf (f) + Sf (0) dt' 

Jo 



for t G [0,5]. Hence, Gronwall's inequality gives 



\(r t ) x \<M5e 



MS 



< M5e 



MS 



in Df. This proves ( |3.2| ) in Df if M\ is sufficiently large and 5 is sufficiently small. 

We next consider the left and right regions Dg, Df which are next to an external end. 
Since the two cases are similar, we will consider the left region only. Let Sj = Si/e for any 
e > 0. Then, the pair (r i; s») satisfies the fixed point equations of either ( [2.1 4j ) or Q2.15Q , 
depending on the type of the boundary condition. Differentiating the equations with respect 
to x and using a slightly modified version of (|3.8| ) where the lower limit a of the integral also 
depends on £: 



d 

d£ 



f{x (t),t) Dg {x {t),t)dt 

= f(x (b) , b) g x (x (b) , b) x^ (b) - f (x (a) , a) g x (x (a) , a) a; € (a) 
— / (a; (a) , a) -Dg (x (a) , a) 



we have 



+ / [/* (x (t) , t) Dg (x (t) ,t)-Df (x (t) , t) g x (x (t) , t)} x c (t) dt, 

(Ci - ifFi - (W) Ut) (0, r) r x + {lf) x U t (x, t) - {lf)U % {xf) x (0, r) 

H jf [«*«). + (W) - (<?). (Wl W).*. (3.10) 

-At) x m,*)+l fitffy.+m) m x -m x (d^)} ( x f) x dt, 

e £ Jo 



where 



Q = 2(u t P l B ) t + E(s 
if the boundary condition is given by (|1.6|) , and 



A 



Of 



R 



A 



if the boundary condition is given by ( |1 . 7| ) . This equation is valid for any e. So, we may 
choose e so small such that 

a=:e|Afr :c (0,t)|max|l, (^4) } < L * e [0, «5] . 



To proceed further, we need an estimate of \r x (0, t) \ . Observe that r (x) satisfies the equation 

xf{r;x,t) = 

where xf (r; x, t) is the solution of the initial value problem 

dx? 



ds 

By differentiation, 



Af (0,r(x))r x + 



\f(xf,s), xf(t;x,t) = x 



dx 

Let Wi = dxf/dx. Then, Wi is the solution to the linear equation 

dw. 



= 0. (3.11) 



ds 

Solving the equation, 



(Af ) (xf (s; x, t) , s) w u Wi(t) = l. 



{Xf) x (xf(s';x,t),s') ds' 



Wi{s) = exp 

Returning to ( |3.11| ), we find 

_1 / r T ( x ) 

t x = —s exp / (Af) (xf (s'\ x, t) , s') ds' 

Af(0,r(x)) y \J t " ; ' J 

Observe that < r (x) < t < 5 and the integrand is bounded. Hence, 

\t x \ < Me MS . (3.12) 

This is the estimate we need. By this estimate, for any Mi, we can choose 5 small enough 
such that the constants a and e are independent of M±. Let Rf (t) and Sf (t) be defined as 
in ( |3.9| ) except that s, is substituted by Sj and D$ is substituted by Dg U Df, We derive 
from (|3.10|) and the identity 

(s t ) t = dfSi - Af (SO, 

that 

Ri (*) + ^ (*) < crSf (t) + M + M5 + M (l + Rf (f) + Sf (0) d* 7 . 

Replacing M and M by M (1 — a) and M (1 — er), respectively, and applying Gronwall's 
inequality, we obtain 

Rf (t) + Sf (t) <[m + Mtf) e MS . 
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Since \sA < \sA, it follows that 



max{\( ri )JA(si)J}<(M + M8) e MS 

in Df U Df. This proves ( |Q| ) in Df U Df if M x is large and 5 is small. 

We next consider the case where the end of the branch is a branching or transitional 
junction. As before, all the branches that are connected to the same junction are considered 
simultaneously. Differentiating the fixed point equation w = Kw where w and Kw are 
defined in ( |2.19| ) and ( |2.2U| ), respectively, we obtain ( |3.1(J| ) in Df U Df for i = j u+ i, . . . , j M 
and 

(h) x = ~ (if). U* (x, t) + \f [(i« Fi ) x + ( W ) (i/o. - (/f), (flf ^)] 

(*), = (*< - tffl - W ) (1, *") ^ + CO. ^ (*, *) - (^), ^ W), (1. T-) ( 3 ' 13 ) 

[(*?#), + (#tf) (^4 - (tf). (flfOi)] W),*, 



in D^ 7 U Df for i = ji, . . . , j v , where 



1=1 F=t/+1 

fl < = e E( m H( 1 ' T ) +e E KAv)t (0 ' r)+ffi 



l'=u+l 



and m* , nj- are defined in the proof of Theorem 2T . Note that the estimate (|3.12j) holds for t x 
in both (|3.10| ) and (|3.13|) , although in the latter case, r is the t-coordinate of the intersection 
of the left-going characteristic curve xf with the vertical line x — 1. The derivation is 
identical. Hence, there is a constant e, independent of M±, such that 



£ N EKWI+ E 
Vfc=i fe'=f+i 

^w(EK(t)l+ E !<,(*) 
\fc=i fe'=i/+i 



< i, 



< i 



in [0, 5}. Let a be the maximum of the quantities on the left hand side of the above inequal- 
ities. Define Rf , Sf 1 , Rf and Sf as in ( |3.9|) with obvious modifications. We see that the 
function 

w (t) = e (t) + s« («)) + e (t) + % (t); 



i=i 



Z'=i/+1 
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satisfies the inequality 

(i--)w(t) < £((i-*)i$(t)+fi£(t))+ £ (i%(*)+(i-^)^(t)) 

i=l l'=v+l 

< M + M5 + M [ (1 + W(t'))dt'. 
Jo 

Hence, by rescaling and using Gronwall's inequality, we achieve 

W(t) < (M + M5^j e MS . 

This proves that 

max{\(r i ) x \,\(s i ) x \}<M 1 
in .Df 1 for i = j 1; . . . , j„ and in Dg for i = j^+i, . . . , j M if Mi is sufficiently large and 5 is 



sufficiently small. We have thus proved (|3.2|) in this case. 

This completes the proof of ( |3.2| ) in all cases. By choosing appropriate values of Mi and 
5, we thus obtain a set Xs t Mo i n the form of ( |3.1|) which is invariant under the mapping S. 

We now show that S is a contraction in Xs,m - Let U = Sv, U = Sv for some v , v £ X5, 
and let W = U — U . W satisfies the vanishing initial and external boundary conditions and 
its differential equations takes the form of (|1.13 ) with the coefficients 



cii = at (x, t, v) , bi = h (x, t, v) , q = q (x, t, v) , A* = A (x, t, v) 
the forcing functions fi and <?, replaced by 

fi =: /j (x, t, v) - fi (x, t, v) + (a* (x, t, u) - a* (x, i, 5)) (3.14) 

and 

dPi 

§i =: # (x, t, v) - gi (x, t, v) + (ft* (x, t, v) - h (x, t, 5)) 



+2 (q (x,t,t>) - Cj (x,t,v)) 



<9x 
3& 



(3.15) 



<9x 

respectively, and the functions Cj in (|2.1|) replaced by 



Cj l — Cj t (x, t, v) — Cj ; (x, t, v) + (Aj ; (x, t, ll) — Aj l (x, t, v)) {Pft — Pjm 

Cj v = Cj v (x, t, v) — Cj v (x, t, v) + [Ajy (x, t, v) — A^, (x, t, t>)) ( Pj Unc — Pj v 
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for I = 1, . . . , u, V = v + 1, . . . , //. By the Lipschitz property and the boundedness 
Mi, there is a constant M such that 



U 



< 



us 



f 



< M\v 



0,5 



1 0,5 ' 



\g\ 0tS <M\v-v\ 0>s , 



C 



< M\v 



0,8 



1 0,5 



Hence, by Theorem 2.1 



\3v - Sv\ o s < MS 



V — V 



0,5 



Therefore, S is contracting in Xs,m if <5 is sufficiently small. 

The rest is standard (cf. e.g., 0). Starting with a v o G X^Mo, we generate an iterative 
sequence Vk+i = Svk- Clearly, each Vk lies in Xs,m and the sequence converges uniformly. 



The limit then satisfies the integral equations in the proof of Theorem [2.1| , and hence, is 
differentiable. Therefore, it is the solution of the quasilinear differential equations. This 
proves the existence and uniqueness of the solution when U 1 = 0. 

If U 1 7^ 0, we regard U 1 as a vector function of x and t and introduce U = U — U 1 . It 
follows that U is a solution of the quasilinear equations (|1.13|) with the forcing functions /$ 
and gi given by 

fi = fi-(Qi) s ai, ~g l = g l -(Pl) x b i -(Q I l ) x 2c l 
and the boundary functions are given by 



P; 



B 



Pf-Pl Qf = Qf -Qi 



and 



C 



Cj, + Pj t , C 



c 



3V 



for I — 1, . . . , v, V — i/+ 1, . . . , fi. Since U has the vanishing initial values, it can be uniquely 
solved for an interval of t e [0, 5]. This gives rise to a solution U. ■ 



Remark: Examples can be constructed to show that if the condition ( |2.4Tj fails at t — 0, 
then, the local solution need not exist or may be not unique. In particular, if fl2.4p fails at 
a source end, then, the system is under-determined, and if it fails at a terminal end, the 
system is over-determined. 

We give next a result for the continuity of dependence of the solution and its derivatives 
on the initial, boundary and forcing functions and their derivatives. This follows from an 



argument similar to the proofs of Lemma |2.1| and Theorem |3J 



Corollary 3.1 Let U = (P,Q) and U = \ P, Qj be two solutions of the quasilinear problem 
of Theorem \3.1\ . Suppose the conditions of that theorem hold for the initial and boundary 
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functions of both solutions. Then, there exists a constant M > 0, independent of initial, 
boundary and forcing functions, such that 



u-u 


< M ( 


pi-pi 


+ 




k,S \ 




C k [0,l] 



Q 1 -Q 1 



+ 5 



f-f 



C k [0,l] 



C k D* 



pB _ pB 



C k [0,S] 



+ 



+ S\g-g 



C k (D s ) 



+ 5 



Q B — Q B 
C-C 



C k [0,8] 



C k [0,S] 

(3.16) 



for k = 0, 1. 



Proof. For k — 0, the result follows from substituting one of the solutions into the coeffi- 
cients, modifying the forcing functions by ( |3.14] )-( |3.15| ), and using Lemma 27T . For k — 1, we 
differentiate the equations and apply the lemma to the resulting equations for the derivatives 
of the solution. The process is standard and is omitted. ■ 
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